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ABSTRACT 

We  present  measurements  of  [Fe/H]  for  six  M  supergiant  stars  and  three  giant  stars  within  2.5  pc  of 
the  Galactic  center  (GC)  and  one  M  supergiant  star  within  30  pc  of  the  GC.  The  results  are  based  on 
high-resolution  (2/A2  =  40,000)  A -band  spectra,  taken  with  CSHELL  at  the  NASA  Infrared  Telescope 
Facility.  We  determine  the  iron  abundance  by  detailed  abundance  analysis,  performed  with  the  spectral 
synthesis  program  MOOG.  The  mean  [Fe/H]  of  the  GC  stars  is  determined  to  be  near  solar,  [Fe/H]  = 
+  0.12  +  0.22.  Our  analysis  is  a  differential  analysis,  as  we  have  observed  and  applied  the  same 
analysis  technique  to  11  cool,  luminous  stars  in  the  solar  neighborhood  with  similar  temperatures  and 
luminosities  as  the  GC  stars.  The  mean  [Fe/H]  of  the  solar  neighborhood  comparison  stars, 
[Fe/H]  =  +0.03  ±  0.16,  is  similar  to  that  of  the  GC  stars.  The  width  of  the  GC  [Fe/H]  distribution  is 
found  to  be  narrower  than  the  width  of  the  [Fe/H]  distribution  of  Baade’s  window  in  the  bulge  but 
consistent  with  the  width  of  the  [Fe/H]  distribution  of  giant  and  supergiant  stars  in  the  solar  neighbor¬ 
hood. 

Subject  headings:  Galaxy:  abundances  —  Galaxy:  center  —  Galaxy:  stellar  content  — 
stars:  abundances  —  stars:  late-type 


1.  INTRODUCTION 

Morris  &  Serabyn  (1996)  and  Serabyn  &  Morris  (1996) 
have  summarized  the  relative  importance  of  star  formation, 
gas  inflow,  and  gas  outflow  in  the  Galactic  center  (GC).  Gas 
from  the  inner  disk  flows  into  the  nucleus,  perhaps  driven 
by  a  bar  potential  (Stark  et  al.  1991;  Morris  &  Serabyn 
1996),  and  some  of  the  stellar  mass  loss  from  the  bulge  may 
fall  into  the  GC  (Blitz  et  al.  1993;  Jenkins  &  Binney  1994). 
Serabyn  &  Morris  (1996)  have  proposed  that  this  inflow  of 
disk  gas  results  in  sustained  star  formation  in  a  “central 
molecular  zone”  (CMZ)  within  a  radius  from  the  GC  of 
R  <  200  pc.  The  GC  is  observed  to  be  currently  forming 
stars  in  the  CMZ  at  a  rate  of  0.3-0.6  Me  yr_1  (Glisten 
1989).  Some  of  the  GC  gas  incorporated  into  stars,  enriched 
by  stellar  nucleosynthesis,  will  be  returned  to  the  GC  inter¬ 
stellar  medium  by  stellar  mass  loss  and  supernovae.  X-ray 
observations  of  hot  gas  in  the  GC  suggest  that  some  gas 
could  be  driven  temporarily  or  permanently  from  the  GC 
by  a  galactic  fountain  or  wind  (Blitz  et  al.  1993;  Morris  & 
Serabyn  1996),  however,  making  it  unclear  how  much 
enriched  material  is  incorporated  into  following  gener¬ 
ations  of  star  formation.  The  presence  of  strong  magnetic 

1  Visiting  Astronomer,  Cerro  Tololo  Inter-American  Observatory. 
CTIO  is  operated  by  AURA,  Inc.  under  contract  to  the  National  Science 
Foundation. 

2  Visiting  Astronomer,  Infrared  Telescope  Facility.  IRTF  is  operated  by 
the  University  of  Hawaii  under  contract  to  the  National  Aeronautics  and 
Space  Administration. 


fields  in  the  GC  also  likely  plays  a  role  in  whether  enriched 
gas  is  driven  from  or  is  retained  within  the  GC. 

Continued  star  formation  in  the  central  few  hundred 
parsecs  of  the  Galaxy  may  lead  to  higher  metallicities 
within  the  CMZ.  Chemical  abundances  in  the  disks  of  spiral 
galaxies  are  observed  to  reach  their  highest  values  at  the 
center  (see  Pagel  &  Edmunds  1981;  Shields  1990).  H  n 
regions,  planetary  nebulae,  and  OB  associations  in  the 
Milky  Way  also  show  a  radial  metallicity  gradient  (Shaver 
et  al.  1983;  Fich  &  Silkey  1991;  Maciel  &  Koppen  1994; 
Simpson  et  al.  1995;  Vilchez  &  Esteban  1996;  Rudolph  et 
al.  1997 ;  Afflerbach,  Churchwell,  &  Werner  1997;  Smartt  & 
Rolleston  1997;  Gummersbach  et  al.  1998).  Chemical  evolu¬ 
tion  models  strive  to  explain  these  gradients  by  considering 
the  relative  star  formation  and  gas  infall/outflow  rates,  and 
the  metal  abundance  of  the  gas  compared  to  the  stars  (see 
Audouze  &  Tinsley  1976;  Rana  1991). 

The  metallicity  of  GC  stars  is  needed  in  order  to  con¬ 
strain  models  and  understand  the  stellar  processes  in  the 
central  parsecs  of  our  Galaxy.  Measurements  of  stellar 
metallicities  in  the  center  of  the  Milky  Way,  which  is 
obscured  by  approximately  30  mag  of  extinction  at  V,  are 
only  now  beginning  to  be  feasible  through  infrared  studies. 
Carr,  Sellgren,  &  Balachandran  (2000)  found  a  solar  Fe 
abundance  in  IRS  7,  an  M  supergiant  at  a  distance  of 
R  =  0.2  pc  from  the  GC,  from  a  detailed  abundance 
analysis  of  CSHELL  (2/A2  =  40,000)  A -band  spectra.  Our 
work  analyzes  similar  data  for  several  M  supergiant  and 
giant  stars  located  in  the  central  cluster,  within  2.5  pc  of  the 
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GC,  and  also  one  M  supergiant  star  located  in  the  Quin¬ 
tuplet  cluster,  30  pc  away  from  the  GC. 

2.  OBSERVATIONS 
2.1.  Sample  Selection 

The  GC  stars  were  selected  with  three  requirements. 
First,  they  have  to  be  brighter  than  K  =  9.5,  a  limit  set  by 
the  sensitivity  of  the  instrument.  An  hour  of  integration 
time  is  required  for  high  signal-to-noise  (S/N)  spectra  for  a 
K  =  9  star  when  observing  with  CSF1ELL  at  the  NASA 
Infrared  Telescope  Facility  (IRTF).  Second,  the  stars  should 
be  cool  luminous  stars  in  order  to  ensure  we  are  observing 
photospheric  absorption  lines  rather  than  the  stellar  wind 
emission  lines  characteristic  of  hot  luminous  stars  in  the 
GC.  Third,  these  cool  luminous  stars  should  have  little  or 
no  water  absorption.  This  is  because  published  line  lists  for 
water  (wavelengths,  excitation  potentials,  dissociation 
potentials,  oscillator  strengths,  and  damping  constants)  are 
not  yet  adequate  for  detailed  abundance  analysis  of  high- 
resolution  stellar  spectra.  The  amount  of  water  absorption 
is  determined  from  H-  and  A -band  low-resolution  spectra 
(X/AX  ~  500)  taken  with  the  IRS  at  CTIO  (S.  V.  Ramirez,  R. 
D.  Blum,  &  K.  Sellgren  2000,  in  preparation).  The  selected 
sample  consists  of  nine  stars  located  in  the  central  cluster 
(R  <  2.5  pc),  and  one  star  located  in  the  Quintuplet  cluster 
(R  =  30  pc).  The  central  cluster  stars  are  IRS  7,  IRS  11,  IRS 
19,  IRS  22,  BSD  72,  BSD  114,  BSD  124,  BSD  129,  and  BSD 
140  (names  and  positions  are  given  by  Blum,  Sellgren,  & 
DePoy  1996a).  The  star  located  in  the  Quintuplet  cluster  is 
VR  5-7  (Moneti,  Glass,  &  Moorwood  1994).  Eleven  cool, 
luminous  stars  in  the  solar  neighborhood  were  selected 
from  the  literature  to  be  observed  and  analyzed  the  same 
way  as  the  GC  stars.  These  11  stars  have  known  abun¬ 
dances  from  detailed  analysis  in  the  optical  (Luck  1982a, 
1982b;  Lambert  et  al.  1984;  Smith  &  Lambert  1985,  1986; 
Luck  &  Bond  1989)  and  in  the  infrared  (Carr  et  al.  2000). 
These  cool,  luminous  stars  in  the  solar  neighborhood  were 


carefully  selected  to  be  in  the  same  range  of  effective  tem¬ 
perature  and  surface  gravity  as  the  GC  stars.  This  is  very 
important  in  order  to  allow  a  differential  analysis  compari¬ 
son,  which  will  cancel  possible  systematic  errors  such  as 
NLTE  effects  (see  §  3.4).  The  selected  cool,  luminous  stars  in 
the  solar  neighborhood  are  listed  in  Table  1. 

2.2.  Data  Acquisition  and  Reduction 

The  R-band  high-resolution  (X/AX  =  40,000)  spectra  used 
in  the  abundance  analysis  were  taken  with  CSHELL  at  the 
IRTF.  The  IRTF  is  a  3  m  telescope  located  at  Mauna  Kea 
Observatory  in  Hawaii,  and  CSHELL  is  the  facility  cryoge¬ 
nic  infrared  echelle  spectrograph  (Tokunaga  et  al.  1990). 
The  observations  were  carried  out  in  1993  May,  1996  June, 
1996  August,  1997  July,  1998  May,  1998  July,  and  1999 
June.  The  detector  for  the  1993  observations  was  a 
256  x  256  NICMOS-3  HgCdTe  array,  which  provided  a 
spectral  coverage  of  « 1000  km  s-1  or  73  A.  The  detector 
for  the  rest  of  the  observations  was  a  256  x  256  SBRC  InSb 
array,  which  gives  a  smaller  coverage  of  x  750  km  s- 1  or  55 
A.  The  slit  used  to  achieve  a  resolution  of  X/AX  =  40,000  is 
0.5"  wide,  30"  long,  with  a  pixel  scale  of  0.25"  pixel-1  for 
1993  data  and  0.20"  pixel-1  for  later  data.  Our  highest 
quality  data  were  taken  in  1998  and  in  1999  when  new 
order-separating  circular  variable  filters  were  installed,  pro¬ 
viding  spectra  free  of  fringes. 

Eight  iron  lines  were  selected  by  inspecting  the  high- 
resolution  atlas  of  cool  stars  (Wallace  &  Hinkle  1996)  and 
by  analyzing  synthetic  spectra,  which  includes  atomic  and 
molecular  features  (see  details  in  §  3  and  in  Carr  et  al.  2000). 
Thousands  of  lines  of  molecules  such  as  CN,  CO  and  H20 
are  present  throughout  the  infrared  spectrum  of  cool  stars. 
The  selected  iron  lines  are  almost  free  of  known  molecular 
contamination  and  unidentified  lines.  The  eight  iron  lines 
are  observed  using  three  grating  settings  of  CSHELL. 
Atomic  parameters  of  the  eight  iron  lines  used  in  the  abun¬ 
dance  analysis  are  listed  in  Table  2. 


TABLE  l 


Stellar  Parameters  for  Nearby  Late-Type  Stars 


Star 

Spectral  Type 

S/Na 

7^  a 
ieff 

(K) 

Mbola 

log  ga 

{“ 

(km  s“*) 

ca 

(km  s-1) 

HR  6146  . 

M6  IIIb 

45° 

3250+  100b 

-5.5” 

0.2  ±0.3” 

2.0  ±0.5' 

9.6  +  1. 1' 

HR  6702  . 

M5  II-IIIb 

36° 

3300+100” 

-3.4” 

0.7  ±0.3” 

2.0  ±0.5' 

7.8  ±1.8' 

HR  7442  . 

M5  IIIasd 

36° 

3450+100“ 

-4.0“ 

0.5  ±0.3“ 

2.4  ±0.2' 

10.5  +  2.7' 

HR  8062  . 

M4  IIIasd 

34° 

3450+100“ 

-3.5“ 

0.7  ±0.3“ 

1.7  ±0.2' 

8.9  ±0.5' 

a  Ori . 

Ml-2  Iab-a' 

305f 

3540  +  260' 

—  7.46 

0.0±0.3g 

2.8  ±0.2' 

14.7  +  0.5' 

HR  8383  . 

M2  Iape  +  h 

100° 

3480  +  250' 

-6.8”“ 

0.0  ±0.3” 

2.7  ±0.2' 

14.4  +  2.1' 

HD  202380 . 

M3  Ibh 

46° 

3600  +  200” 

1 

p* 

0.6  ±0.5” 

2.5  ±0.2' 

16.2  +  0.5' 

HD  163428 . 

K5  IIb 

67° 

3800  +  200” 

-5.5”“ 

0.6  ±0.5” 

2.3  ±0.2' 

11.5  +  3.2' 

BD  +  59°594 . . . . 

Ml  Ibj 

49° 

4000  ±  200j 

-6.6* 

—  0.9  ±0.3* 

3.1  ±0.2' 

18.2+1.0' 

HD  232766 . 

Ml  lab* 

38° 

4000  ±  200j 

-6.6* 

0.2  ±0.3* 

2.2  ±0.2' 

11.5  +  3.7' 

HR  8726  . 

K5  Iabk 

57° 

4000±200k 

—  5.2k* 

0.5±0.3k 

2.4  ±0.2' 

11.4  +  0.5' 

a  S/N  =  mean 

signal-to-noise  ratio 

per 

pixel;  Teff  = 

effective  temperature;  Mbol 

=  bolometric  magnitude; 

g  =  surface  gravity;  {  =  microturbulent  velocity;  f  =  macro  turbulent  velocity. 
”  Smith  &  Lambert  1985. 

'  This  paper. 

“  Smith  &  Lambert  1986. 

'  Carr  et  al.  2000. 
f  Wallace  &  Hinkle  1996. 

B  Lambert  et  al.  1984. 
h  Luck  1982b. 

1  Lang  1991,  p.  142. 
j  Luck  &  Bond  1989. 
k  Luck  1982a. 
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TABLE  2 
Fe  Line  Data 


Wavelength 

(A) 

X 

(eVf 

log  gf* 

21781.82 . 

3.415 

-4.485 

22381.27 . 

5.844 

-1.458 

22386.90 . 

5.033 

-0.481 

22391.22 . 

5.320 

- 1.600 

22398.98 . 

5.099 

- 1.249 

22818.82 . 

5.792 

-1.296 

22838.60 . 

5.099 

-1.325 

22852.17 . 

5.828 

-0.612 

a  x  =  excitation  potential;  gf  =  gf- 
value  determined  in  this  paper. 


Spectra  were  acquired  at  two  different  positions  along  the 
slit  and  reduced  separately  in  order  to  determine  any  sys¬ 
tematic  effects  such  as  fringes.  Several  individual  spectrum 
pairs  were  observed  per  star,  until  the  desired  signal-to- 
noise  ratio  was  reached.  Stars  of  spectral  type  A  or  B  were 
observed  as  close  to  the  program  star’s  air  mass  as  possible 
in  order  to  correct  for  telluric  absorption  features  and  to 
remove  fringes.  Such  stars  have  no  significant  spectral  fea¬ 
tures  in  the  observed  wavelength  regions. 

Image  Reduction  and  Analysis  Facility  (IRAF3)  and 
VISTA  were  used  for  data  reduction.  The  reduction  process 
began  by  flat-fielding  the  individual  spectra  with  calibration 
lamp  flats,  which  were  taken  for  every  grating  setting.  Sky 
subtraction  was  performed  by  subtracting  two  consecutive 
spectra  taken  at  different  positions  along  the  slit.  Bad  pixels 
were  replaced  by  an  interpolated  value  computed  from 
neighboring  pixels  in  both  the  dispersion  and  perpendicular 
directions.  Individual  spectra  were  extracted  with  a  5  pixel 
wide  aperture  using  the  APSUM  package  in  IRAF. 
Extracted  spectra  observed  at  the  same  position  along  the 
slit  were  averaged  to  produce  two  spectra  per  star,  one 
spectrum  for  each  position  along  the  slit.  Spectra  of  the 
program  stars  were  then  divided  by  the  appropriate  A-  or 
B-type  atmospheric  standards,  observed  at  the  same  posi¬ 
tion  along  the  slit  and  reduced  in  the  same  way,  to  remove 
telluric  absorption  features  and  fringes.  Each  spectrum  was 

3  The  IRAF  software  is  distributed  by  the  National  Optical  Astronomy 
Observatories  under  contract  with  the  National  Science  Foundation. 


then  multiplied  by  a  blackbody  of  the  same  temperature  as 
the  A  or  B  star  to  put  the  spectra  on  a  relative  flux  density 
(Fj  scale.  The  temperature  of  the  A  or  B  star  was  deter¬ 
mined  from  its  spectral  type.  More  than  five  telluric  absorp¬ 
tion  lines  were  used  to  obtain  wavelength  solutions.  The 
wavelength  calibration  was  performed  using  the  tasks  of 
VISTA.  Spectra  were  shifted  in  wavelength  to  correct  for 
radial  velocity  differences,  by  comparing  the  observed  and 
synthetic  spectra.  The  spectra  were  normalized  by  a  linear 
fit  of  continuum  bands.  These  continuum  bands  were 
chosen  by  comparing  the  observed  spectrum  of  a  Ori  from 
the  high-resolution  atlas  (Wallace  &  Hinkle  1996)  with  the 
synthetic  spectrum  of  a  Ori.  The  bands  selected  are  regions 
in  the  spectrum  in  which  the  synthesis  and  the  normalized 
observations  of  a  Ori  are  coincident  and  equal  to  unity.  We 
estimated  the  error  per  pixel  as  the  difference  of  obser¬ 
vations  taken  at  two  positions  along  the  slit.  We  computed 
the  mean  signal-to-noise  ratio  (S/N)  per  pixel  by  the  average 
of  the  signal-to-noise  ratio  per  pixel  for  all  spectra  observed 
for  that  star.  The  mean  S/N  per  pixel  for  each  star  is  listed  in 
Tables  1  and  3. 

3.  ABUNDANCE  ANALYSIS 

The  abundance  analysis  was  done  using  a  current  version 
of  the  LTE  spectral  synthesis  program  MOOG  (Sneden 
1973).  The  program  requires  a  line  list  giving  the  wave¬ 
length,  excitation  potential,  (//-values,  and  damping  con¬ 
stants  for  all  atomic  and  molecular  lines  that  contribute  to 
the  spectrum.  The  method  used  to  construct  the  line  list  is 
described  in  §  3.1.  In  addition,  an  input  model  atmosphere 
for  the  effective  temperature  and  surface  gravity  appropri¬ 
ate  for  each  star  and  a  value  for  the  microturbulent  velocity 
is  also  required.  The  solar  abundance  model  atmospheres 
from  Plez  (Plez  1992;  Plez,  Brett,  &  Nordlund  1992)  were 
used  in  our  analysis  because  they  are  the  only  ones  that 
include  gravities  and  temperatures  as  low  as  those  of  the 
Galactic  center  stars.  The  Plez  models  also  include  spher¬ 
icity,  which  is  appropriate  for  supergiant  stars  (Plez  et  al. 
1992).  MOOG,  however,  does  not  account  for  sphericity  in 
its  calculations.  Carr  et  al.  (2000),  who  also  used  MOOG, 
compute  the  effective  temperature  and  microturbulent 
velocity  of  a  Ori,  VV  Cep  (HR  8383),  and  [l  And  using 
different  sets  of  model  atmospheres,  including  the  Plez 
models  and  two  plane-parallel  atmospheres.  When  the 


TABLE  3 

Stellar  Parameters  for  Galactic  Center  Stars 


Ma  p  c 

Star  Luminosity  Class  S/Ns  (K)  Mbols,b  (M0)  log  g*  (km  s  3)  (km  s  *) 


IRS  7°  .  I  52  3470  +  250  -9.0  17  +  3  -0.6  ±0.2  3.3  ±0.4  20.6  +  2.7 

VR  5-7 .  I  30  3500  +  300  -7.8  14  +  2  -0.2 ±0.3  2.9  +  0.5  12.6+1.6 

IRS  19 .  I  75  3650  +  300  -7.2  14  +  2  0.1  ±0.3  2.7  +  0.5  13.4  +  2.0 

IRS  22 .  I  25  3550  +  300  -6.4  10  +  2  0.3  +  0.3  2.5  +  0.5  12.8  +  1.6 

BSD  124 .  I  10  3600  ±300  -5.5  7  +  3  0.4  ±0.3  2.4  ±0.5  12.7  +  2.7 

BSD  129 .  I  13  3650  +  300  -5.3  7  +  3  0.5  +  0.3  2.4  +  0.5  12.3  +  2.7 

BSD  72 .  I  13  3750  +  300  -4.5  5  +  2  0.8  +  0.3  2.1  +  0.5  11.0  +  2.7 

BSD  114 .  Ill  15  3100  +  280  -5.8  3  +  1  -0.2  ±0.5  2.9  +  0.6  9.0  +  3.0 

IRS  11 .  Ill  16  3100  +  280  -5.3  3  +  1  0.0  +  0.5  2.8  +  0.6  9.1  ±2.2 

BSD  140 .  Ill  13  3100  +  280  -4.8  2.5+1  -0.1  ±0.5  2.9  +  0.6  9.0  +  3.0 


a  S/N  =  mean  signal-to-noise  ratio  per  pixel  ;  Gr  =  effective  temperature;  Mbol  =  bolometric  magnitude;  M  =  mass; 
g  =  surface  gravity;  {  =  microturbulent  velocity;  f  =  macroturbulent  velocity. 

b  Error  in  Mbol  is  ±  0.4  (Carr  et  al.  2000),  dominated  by  the  uncertainty  in  the  extinction  curve  \E(H  —  K)/A  J. 
c  Stellar  parameters  from  Carr  et  al.  2000. 
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stellar  parameters  derived  from  each  model  were  used  self- 
consistently  to  derive  [Fe/H],  the  result  for  [Fe/H]  was  the 
same.  The  determination  of  the  stellar  parameters  for  our 
sample  is  discussed  in  §  3.2. 

3.1.  Atomic  and  Molecular  Parameters 

The  line  list  was  created  the  same  way  as  that  of  Carr  et 
al.  (2000),  using  two  main  steps:  (1)  The  solar  spectrum  was 
used  to  determine  the  (//-values  and  the  damping  constants, 
and  (2)  minor  adjustments  were  made  by  comparison  to  the 
Arcturus  spectrum  (Wallace  &  Hinkle  1996).  An  initial  line 
list  was  compiled  using  Fe  I  line  positions  and  energy  levels 
from  Nave  et  al.  (1994)  and  the  Kurucz  (1993)  line  lists  for 
CN  and  other  atomic  lines.  The  damping  constants  for  all 
atomic  lines  were  initially  set  to  twice  that  of  the  Unsold 
approximation  for  van  der  Waals  broadening  (Holweger  et 
al.  1991).  A  synthetic  spectrum  was  generated  for  the  Sun, 
using  the  Kurucz  (1993)  solar  model.  The  .(//-values  and 
damping  constants  were  adjusted  to  match  the  observed 
solar  spectrum  (Livingstone  &  Wallace  1991),  when  needed, 
and  unidentified  lines  were  noted.  Then,  a  synthetic  spec¬ 
trum  was  generated  for  Arcturus,  using  the  model  atmo¬ 
sphere,  stellar  parameters,  and  abundances  from  Peterson, 
Dalle  Ore,  &  Kurucz  (1993).  The  solar  line  list  provided  a 
good  match  to  the  Arcturus  spectrum.  The  gf- values  of  low- 
excitation  lines,  not  observable  in  the  Sun,  had  to  be  adjust¬ 
ed.  The  final  Fe  i  atomic  parameters  are  listed  in  Table  2.  A 
synthetic  spectrum  of  a  Ori,  generated  using  the  final  line 
list,  is  compared  in  Figure  1  to  the  observed  spectrum  from 
Wallace  &  Hinkle  (1996)  in  the  region  of  each  of  our  eight 


Fe  i  lines.  A  Plez  model  atmosphere  with  the  stellar  param¬ 
eters  listed  in  Table  3,  and  the  CNO  abundances  from 
Lambert  et  al.  (1984)  were  used. 

The  depth  of  line  formation  was  examined  using  contri¬ 
bution  functions,  which  provide  an  indication  of  spectral 
line  formation  coming  from  different  layers  of  a  stellar 
model  atmosphere  (see  details  in  Edmonds  1969;  Sneden 
1973).  The  determination  of  the  depth  of  line  formation  is 
needed  to  make  sure  each  line  is  formed  in  the  range  of 
opacity  (t)  covered  by  our  model  atmosphere.  The  eight 
iron  lines  used  in  this  abundance  analysis  are  formed  in  the 
range  of  —  2.0  <  log  t  <  0.3,  where  Plez  model  atmo¬ 
spheres  cover  a  range  of  —  5.6  <  log  r  <  2.6.  Note  that  the 
Sc  i  lines  in  supergiant  stars,  which  neither  we  nor  Carr  et 
al.  (2000),  who  included  hyperline  splitting  calculations  for 
Sc  i  lines,  are  able  to  model  correctly,  are  shown  by  the 
contribution  function  to  be  partially  formed  at  optical 
depths  outside  the  range  of  the  Plez  models.  We  therefore 
approximately  model  Sc  i  lines  by  scaling  the  Fe  i  line  pro¬ 
files  to  the  correct  depth  and  width  to  allow  us  to  separate 
the  one  Fe  i  line  at  22399  A  that  is  blended  with  a  Sc  i  line  at 

o 

22400  A.  In  the  case  of  IRS  7  this  procedure  was  not  suc¬ 
cessful. 

3.2.  Stellar  Parameters 
3.2.1.  Effective  Temperature 

The  effective  temperature  (Teff)  is  a  key  parameter  for  any 
abundance  analysis.  For  the  cool,  luminous  stars  in  the 
solar  neighborhood,  the  effective  temperatures  were  taken 
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Fig.  1. — Observed  spectrum  (filled  squares )  of  a  Ori  from  Wallace  &  Hinkle  (1996)  compared  with  synthetic  spectra  generated  by  MOOG,  using  CNO 
abundances  from  Lambert  et  al.  (1984),  model  atmosphere  from  Plez  (1992),  and  stellar  parameters  listed  in  Table  1.  The  error  bars  come  from  the  difference 
of  observations  taken  in  two  epochs.  We  overplot  a  synthetic  spectrum  (thick  line )  derived  using  the  Fe  abundances  listed  in  Table  5  and  synthetic  spectra 
(thin  lines )  computed  with  Fe  abundances  different  by  +0.20  dex  ,  which  corresponds  to  the  typical  uncertainty  in  [Fe/H]  for  a  solar  neighborhood  cool, 
luminous  star.  Fe  lines  and  their  wavelengths  in  A  are  marked  ( bold  vertical  lines )  at  the  top  of  each  panel.  CN,  Sc,  and  unidentified  lines  (question  marks )  are 
also  marked  (vertical  lines).  Tick  marks  along  the  x-axis  are  1  A  apart.  The  [Fe/H]  value  from  Table  5  for  each  line  is  given  in  each  panel. 
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from  the  literature.  Some  effective  temperatures  come  from 
spectroscopic  measurements  (Luck  1982a,  1982b;  Luck  & 
Bond  1989;  Carr  et  al.  2000).  Other  effective  temperatures 
come  from  the  calibration  of  V  —  K  colors  versus  Tett,  where 
T*t  comes  from  angular  diameters  measured  by  infrared 
interferometry  (Smith  &  Lambert  1985,  1986).  a  Ori  has  a 
Te{{  determination  from  its  angular  diameter  measured  by 
infrared  interferometry  (Tcff  =  3605  +  43  K;  Dyck,  van 
Belle,  &  Thompson  1998),  which  is  very  good  agreement 
with  the  spectroscopic  value  obtained  by  Carr  et  al.  (2000) 
(Teff  =  3540  ±  260;  see  Table  1). 

For  IRS  7,  the  brightest  infrared  source  in  the  GC,  Carr 
et  al.  (2000)  determined  Tcff  spectroscopically  by  requiring 
the  carbon  abundance  derived  from  CO  lines  to  be  indepen¬ 
dent  of  the  lower  excitation  potential.  We  adopt  the  stellar 
parameters  for  IRS  7  from  Carr  et  al.  (2000),  which  were 
derived  using  the  standard  Plez  grid  of  model  atmospheres. 
Their  Teff  analysis  relied  largely  on  the  second-overtone  CO 
lines  present  in  the  77  band  (1.5-1. 8  lira),  because  most  of  the 
first-overtone  CO  lines  present  in  the  K  band  are  saturated. 

For  the  remaining  nine  stars  in  the  GC,  the  77-band 
spectra  are  too  faint  to  be  observed  with  CSHELL  at  the 
IRTF  because  of  high  extinction.  Therefore  CO  cannot  be 
used  to  compute  the  effective  temperature.  For  those  nine 
GC  stars,  broad  molecular  bands  present  in  77-  and  TC-band 
low-resolution  spectra  (2/A2  ~  500)  are  used  to  estimate 
Teff  (S.  V.  Ramirez,  R.  D.  Blum,  &  K.  Sellgren  2000,  in 
preparation).  These  features  are  CO  (2.3  /im)  and  H20  (1.9 
yum).  The  CO  strength  increases  with  decreasing  Tett  and 
decreasing  surface  gravity  (log  g).  The  H20  strength  also 
increases  with  decreasing  Fcff  but  increases  with  increasing 
log  g.  These  two  features  together,  therefore,  provide  two- 
dimensional  spectral  classification  (Baldwin,  Frogel,  & 
Persson  1973;  Kleinmann  &  Hall  1986;  Blum,  Sellgren  & 
DePoy  1996b).  The  low-resolution  spectra  of  GC  stars  were 
taken  with  the  IRS  at  CTIO,  used  in  cross-dispersed  mode 
to  acquire  the  77  and  K  bands  simultaneously  (S.  V. 
Ramirez,  R.  D.  Blum,  &  K.  Sellgren  2000,  in  preparation). 
Nearby  late-type  stars  with  Tm  determinations  based 
directly  or  indirectly  on  the  infrared  interferometry  tech¬ 
nique  (Smith  &  Lambert  1985,  1986,  1990;  Fernandez- 
Villacanas,  Rego,  &  Cornide  1990;  McWilliam  1990;  Dyck 
et  al.  1998;  Richichi  et  al.  1998)  were  observed  to  calibrate 
our  determination.  The  low-resolution  spectra  of  cool, 
luminous  stars  with  known  were  taken  with  the  IRS  at 
CTIO  and  with  TIFKAM  at  MDM  Observatory  (S.  V. 
Ramirez,  R.  D.  Blum,  &  K.  Sellgren  2000,  in  preparation). 
The  distinction  between  giants  and  supergiants  is  deter¬ 
mined  by  the  presence  or  weakness  of  H20,  respectively. 
The  CO  index  is  computed  in  the  same  way  as  in  Blum  et  al. 
(1996b).  Once  the  distinction  between  giants  and  super¬ 
giants  is  established,  is  determined  from  the  relations 
shown  in  Figure  2,  which  are  the  best  unweighted  linear  fit 
to  the  data.  VR  5-7,  IRS  19,  IRS  22,  BSD  72,  BSD  124,  and 
BSD  129  were  classified  as  supergiant  stars,  and  IRS  11, 
BSD  114,  and  BSD  140  were  classified  as  giant  stars.  The 
luminosity  class  for  GC  stars  is  also  listed  in  Table  3.  Our 
values  for  Fcff  for  GC  stars  agree  with  the  results  from  Blum 
et  al.  (1996b),  for  the  stars  that  are  common  to  both 
samples.  Our  value  for  reff  from  the  CO  index  of  IRS  7, 
3400  ±  300  K,  is  also  in  good  agreement  with  the  spectro¬ 
scopic  value  of  Teff  from  Carr  et  al.  (2000),  3470  ±  250  K. 
IRS  7  has  a  low  carbon  abundance  (Carr  et  al.  2000,  [C/ 
H]  =  —  0.8).  The  good  agreement  between  the  Teff  from  CO 


Fig.  2. — CO  index  vs.  effective  temperature,  Teff,  for  supergiant  and 
giant  stars.  Upper  panel :  The  calibration  for  supergiant  stars,  with  T* 
from  Dyck  et  al.  (1998)  (filled  squares )  and  Richichi  et  al.  (1998)  (open 
circles),  both  measured  from  infrared  interferometry.  Lower  panel :  The 
calibration  for  giant  stars.  Values  of  Tetf  from  McWilliam  (1990)  (open 
triangles),  Smith  &  Lambert  (1985,  1986,  1990)  (filled  circles),  and 
Fernandez-Villacanas  et  al.  (1990)  (open  stars)  are  derived  from  a  relation¬ 
ship  between  V  —  K  and  ^eff  measured  by  infrared  interferometry.  Teff  from 
Dyck  et  al.  (1998)  (open  circles)  are  derived  from  infrared  interferometry 
measurements.  In  both  panels  the  uncertainties  come  from  the  literature. 
The  best  unweighted  linear  fit  to  the  data  (solid  lines)  is  used  to  compute 
the  Ttff  of  the  Galactic  center  stars. 

index  and  the  spectroscopic  Teff  suggests  that  our  derived 
values  of  Te  ff  are  not  strongly  affected  by  moderate  varia¬ 
tions  in  the  abundance  of  carbon.  The  typical  uncertainty  in 
Teff  computed  from  the  CO  index  is  300  K  for  supergiants 
and  280  K  for  giants.  The  Teff  for  cool,  luminous  stars  in  the 
solar  neighborhood  is  listed  in  Table  1  and  for  GC  stars  is 
listed  in  Table  3. 

3.2.2.  Surface  Gravity 

For  the  cool,  luminous  stars  in  the  solar  neighborhood, 
the  surface  gravity,  g,  from  the  literature  was  used  (Luck 
1982a,  1982b;  Smith  &  Lambert  1985,  1986;  Luck  &  Bond 
1989;  Carr  et  al.  2000).  For  the  GC  stars,  the  surface  gravity 
is  determined  from  the  relation 

log  g  =  log  (M/Mq)  +  4  log  (Teff/Teff0) 

-  log  (L/Lq)  -  log  gQ  .  (1) 

T*  f  has  been  already  determined.  The  luminosity,  L,  is  com¬ 
puted  from  Mbol,  which  is  determined  from  K,  AK,  the  GC 
distance,  and  the  bolometric  correction,  BC.  A  GC  distance 
of  8  kpc  is  assumed  (Reid  1993),  K  and  AK  are  from  Blum  et 
al.  (1996a),  and  the  BC  is  taken  from  Elias,  Frogel,  & 
Humphreys  (1985).  The  final  Mbol  for  the  GC  stars  is  listed 
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in  Table  3.  The  uncertainty  in  Mbol  from  the  GC  stars 
comes  mainly  from  uncertainties  in  the  extinction  curve, 
E(H —K)/Ak  ,  and  it  is  estimated  to  be  ±0.4  (Carr  et  al. 
2000).  Once  Mbol  and  Teff  are  determined,  the  stars  are 
placed  in  the  H-R  diagram,  and  masses,  M,  are  obtained  by 
overplotting  evolutionary  tracks  of  varying  initial  masses. 
Solar  metallicity  evolutionary  tracks  (Schaller  et  al.  1992) 
are  used  to  determine  masses  for  supergiant  stars.  For  the 
coolest  stars,  that  fall  outside  the  Tm  range  of  Schaller  et  al. 
(1992)  evolutionary  tracks,  asymptotic  giant  branch  (AGB) 
tracks  from  Marigo,  Bressan,  &  Chiosi  (1996)  were  used  to 
estimate  their  masses.  Note  that  the  stellar  evolution 
models  assume  mass  loss  for  the  more  massive  stars,  and  the 
current  mass  of  the  star  at  its  present  age  is  assumed.  Twice 


o 


uo 

o 


solar  metallicity  evolutionary  tracks  (Schaerer  et  al.  1993) 
gave  the  same  result  for  the  derived  masses.  The  value  of  log 
g  for  GC  stars  is  listed  in  Table  3.  The  uncertainty  is  esti¬ 
mated  considering  the  uncertainties  in  mass,  effective  tem¬ 
perature,  and  luminosity.  Figure  3  shows  the  FI-R  diagram 
of  GC  stars  and  also  the  cool,  luminous  stars  in  the  solar 
neighborhood.  It  is  seen  that  the  GC  stars  and  the  cool, 
luminous  stars  in  the  solar  neighborhood  occupy  the  same 
place  in  the  H-R  diagram  and  hence  are  very  similar  types 
of  stars. 

3.2.3.  Microturbulent  Velocity 

The  values  for  the  microturbulent  velocity  (£)  for  the  cool, 
luminous  stars  found  in  the  literature  show  considerable 


lQg  (Teff) 

Fig.  3. — H-R  diagram  of  Galactic  center  (GC)  stars  (filled  circles)  and  cool,  luminous  stars  in  the  solar  neighborhood  ( open  squares ).  Stellar  parameters 
and  their  uncertainties  are  given  in  Tables  1  and  3.  Solar  metallicity  evolutionary  tracks  (Schaller  et  al.  1992)  are  overplotted  to  determine  masses  for  GC 
stars.  The  GC  stars  and  the  cool,  luminous  stars  in  the  solar  neighborhood  occupy  the  same  place  in  the  H-R  diagram  and  hence  are  very  similar  types  of 
stars. 
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variation  among  authors.  Carr  et  al.  (2000)  found  for  VV 
Cep  (HR  8383)  a  much  lower  value  ({  =  3.7  ±  0.2  km  s-1) 
than  the  one  from  Luck  (1982b),  £  =  5.0  ±  0.5  km  s” 1.  Carr 
et  al.  (2000)  used  CO  lines  of  the  same  excitation  potential, 
for  which  the  carbon  abundance  should  be  independent  of 
the  equivalent  widths  of  the  lines  when  £  is  correct.  Luck 
(1982b)  used  the  same  principle  but  used  Fe  i  lines  instead  of 
CO  lines. 

Seven  of  our  Fe  i  lines  have  similar  excitation  potential, 
and  those  lines  were  used  to  derive  the  microturbulence  for 
the  cool,  luminous  stars  in  the  solar  neighborhood.  This 
was  done  in  an  iterative  way.  An  Fe  abundance  was  derived 
from  the  synthetic  model  of  cool,  luminous  stars  in  the  solar 
neighborhood,  which  included  lines  of  Fe,  CN,  and  other 
lines.  Using  the  derived  Fe  abundance,  a  synthetic  model 
was  computed  containing  only  Fe  lines,  and  an  equivalent 
width  was  measured  for  the  Fe  lines.  This  way  the  equiva¬ 
lent  width  is  free  from  contributions  from  CN  or  other 
species.  Then,  the  derived  Fe  abundances  from  the  synthesis 
were  plotted  versus  the  Fe  equivalent  widths.  This  process 
was  repeated  with  different  microturbulent  velocities  until 
the  Fe  abundance  was  independent  of  the  Fe  equivalent 
width.  The  final  equivalent  widths  of  the  Fe  lines  used  in 
this  process  are  listed  in  Table  4.  The  obtained  value  of  £  for 
cool,  luminous  stars  in  the  solar  neighborhood  is  listed  in 
Table  1.  The  uncertainty  in  the  obtained  microturbulent 
velocity  comes  from  the  uncertainty  in  the  slope  in  the 
graph  of  Fe  abundance  versus  equivalent  width,  based  on 
the  scatter  of  the  data  points.  Our  value  of  £  for  a  Ori 
(£  =  2.8  ±  0.2)  is  slightly  different  from  the  value  of  £  from 
Carr  et  al.  (2000)  (£  =  3.23  ±  0.15).  There  is  good  agreement 
between  our  values  of  £  and  the  ones  given  by  Smith  & 
Lambert  (1985,  1986),  where  the  mean  difference  is  0.15  km 
s-1.  Our  values  of  £  are  systematically  lower  than  the 
values  given  by  Luck  (1982b)  and  Luck  &  Bond  (1989)  (HR 
8383,  HD  202380,  HD  163428,  BD  +  59°594,  HD  232766), 
where  the  mean  difference  is  2.4  km  s-1.  Our  value  of  £  for 
HR  8726  is  in  good  agreement  with  the  value  given  by  Luck 
(1982a). 

We  do  not  have  sufficient  S/N  for  the  GC  stars  to  apply 
this  technique  for  finding  £.  Instead,  to  get  the  micro¬ 
turbulent  velocity  for  the  GC  stars,  a  relation  between  log  g 
and  £  is  used.  McWilliam  (1990)  showed  that  log  g  and  £  are 
related  in  G  and  K  giants.  The  values  of  £  and  log  g  for  cool, 
luminous  stars  in  the  solar  neighborhood,  including  the 


log  g 

Fig.  4. — Relation  between  microturbulent  velocity  (£)  and  surface 
gravity  (log  g).  Values  of  £  and  log  g  are  from  Table  1  of  this  paper  (open 
squares),  Carr  et  al.  (2000)  (crosses),  and  Smith  &  Lambert  (1985,  1986, 
1990)  (filled  triangles).  The  best  unweighted  linear  fit  to  the  data  ( solid  line) 
is  used  to  compute  {  for  the  GC  stars. 

values  from  Carr  et  al.  (2000)  and  Smith  &  Lambert  (1985, 
1986,  and  1990),  are  plotted  to  get  a  relation  to  apply  to  the 
GC  stars.  Figure  4  shows  this  relation  for  stars  with  Te{{ 
and  log  g  similar  to  those  of  the  GC  stars.  An  unweighted 
linear  fit  to  the  data  gives 

(  =  (2.78  -  0.82  x  log  g)  km  s”1  .  (2) 

The  uncertainty  in  the  microturbulent  velocity  is  esti¬ 
mated  from  the  uncertainty  of  the  fit  ( ±  0.4  km  s  " 1)  and  the 
uncertainty  from  the  surface  gravity.  The  value  of  £  for  IRS 
7  obtained  by  this  fit,  {  =  3.3  ±  0.4  km  s-1,  is  in  good 
agreement  with  the  value  of  £  from  Carr  et  al.  (2000), 
£  =  3.30  ±  0.34  km  s_  1,  obtained  using  CO  lines. 

3.2.4.  Macroturbulent  Velocity 

The  synthetic  spectra  have  to  be  convolved  with  a  macro- 
turbulent  broadening  function  (Gray  1992,  p.  405),  and  with 
instrumental  broadening,  to  match  the  line  widths  of  the 


TABLE  4 

Fe  n  Equivalent  Widths  for  Nearby  Late-Type  Stars8 


Wavelength  of  Fe  i  Lines  (A) 


Star 

22381.3 

22386.9 

22391.2 

22399.0 

22818.8 

22838.6 

22852.2 

HR  6146  . 

90  +  6 

355  +  8 

162  +  7 

210  +  7 

130  +  6 

277  +  7 

211  +  7 

HR  6702  . 

67  +  5 

332  +  9 

136  +  7 

189  +  7 

119  +  6 

234  +  7 

220  +  7 

HR  7442  . 

60  +  5 

339  +  9 

143  +  8 

174  +  8 

141  +  8 

294  +  8 

246  +  8 

HR  8062  . 

65  +  5 

351  +  9 

149  +  7 

199  +  6 

110  +  6 

208  +  7 

a  Ori . 

101  +  8 

520  +  10 

218+10 

319  +  10 

178  +  10 

339  +  10 

271  ±  10 

HD  202380 . 

102  +  8 

476  +  11 

180  +  9 

311  +  9 

137  +  8 

299  +  9 

HR  8383  . 

86  +  7 

466  +  10 

183  +  9 

285  +  9 

112  +  8 

261  +  9 

233  +  9 

HD  163428 . 

87  +  7 

430  +  9 

172  +  9 

241  +  9 

113  +  8 

263  +  8 

BD  +59°594 . 

107  +  9 

522  +  11 

171  +  12 

329  +  12 

96  +  9 

289  +  13 

226  +  12 

HD  232766 . 

66  +  6 

410  +  8 

159  +  9 

228  +  9 

97  +  8 

216  +  9 

173  +  9 

HR  8726  . 

71  +  6 

397  +  8 

158  +  8 

234  +  8 

115  +  8 

239  +  8 

186  +  8 

a  Equivalent  widths  of  each  Fe  I  line  (mA). 
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observed  spectrum.  The  instrumental  broadening  is  a 
Gaussian  of  full  width  half-maximum  given  by  the  spectral 
resolution  of  the  instrument  (7  km  s”1).  CN  lines  in  two 
bands  (21798-21804  A;  22403-22409  A)  were  used  to  deter¬ 
mine  the  macroturbulent  velocity  (£)  for  the  cool,  luminous 
stars  in  the  solar  neighborhood.  Several  synthetic  spectra 
were  generated  by  MOOG  using  different  values  of  £  and  N 
abundances.  Then,  x2  was  computed  for  each  synthetic 
spectrum,  and  the  best  set  of  N  abundance  and  £  was 
chosen  for  the  minimum  x2-  Note  that  this  procedure  is  not 
measuring  £  but  rather  is  determining  the  best  line  broaden¬ 
ing  that  fits  the  data.  Also,  the  N  abundance  obtained  by 
this  technique  is  not  real  because  the  C  and  O  abundances 
have  not  been  measured  independently.  Figure  5  shows  the 
X2  contour  maps,  one  for  each  band,  for  HR  8726.  Figure  5 
also  shows  the  observed  spectrum  of  HR  8726  and  the  syn¬ 
thetic  spectrum  for  the  values  of  £  and  N  abundance  that 
give  a  minimum  x2  for  both  bands.  The  macroturbulent 
velocity  for  the  cool,  luminous  stars  in  the  solar  neighbor¬ 
hood  is  given  in  Table  1.  The  uncertainty  in  £  is  given  by  the 
standard  deviation  of  the  values  of  £  derived  for  each  band. 


This  uncertainty  is  more  conservative  than  the  uncertainty 
derived  from  x2  statistics. 

For  the  GC  stars,  the  same  method  was  used  for  IRS  7, 
VR  5-7,  IRS  19,  IRS  22,  and  IRS  11.  For  the  remaining 
stars,  BSD  72,  BSD  114,  BSD  124,  BSD  129,  and  BSD  140, 
only  one  band  with  CN  lines  was  observed,  and  the  band 
did  not  have  enough  S/N  to  carry  out  this  method.  The 
values  of  £  in  Table  1  were  used  to  derive  a  relationship 
between  log  g  and  £  for  supergiant  stars.  Gray  &  Martin 
(1979)  established  a  relationship  between  £  and  £  for  giant 
stars.  McWilliam  (1990)  shows  that  £  depends  on  log  g,  so 
by  combining  these  two  dependences  we  derive  a  relation 
between  £  and  log  g.  The  values  of  £  for  supergiant  stars 
measured  by  the  broadening  of  the  CN  lines  are  plotted 
versus  log  g  to  get  a  relation  to  apply  to  the  fainter  GC 
supergiant  stars.  Figure  6  shows  this  relation.  An 
unweighted  linear  fit  to  the  data  gives 

£  =  (14.5  —  4.40  x  log  g)  km  s_1  .  (3) 

The  uncertainty  in  the  microturbulent  velocity  is  estimated 
from  the  uncertainty  of  the  fit  ( ±  2.4  km  s“ 1)  and  the  uncer- 


£(km/ s) 


Mk) 


Fig.  5. — Left-hand  panels:  y2  contour  maps  for  determining  the  nitrogen  abundance,  log  t  (N),  and  macroturbulent  velocity  (()  for  HR  8726.  The 
minimum  y2  gives  the  best  set  of  log  e  (N)  and  f  for  each  set  of  CN  lines  (21798-21804  A  at  the  top,  and  22403-22409  A  at  the  bottom).  Right-hand  panels: 
Observed  spectrum  of  HR  8726  (filled  squares)  with  error  bars  from  the  difference  of  observations  taken  at  two  positions  along  the^slit  and  the  synthetic 
spectrum  of  HR  8726  (solid  line)  using  the  derived  log  €  (N)  and  £  for  each  CN  band  (21798-21804  A  at  the  top,  and  22403-22409  A  at  the  bottom).  The 
synthetic  spectrum  is  calculated  by  MOOG,  using  stellar  parameters  in  Table  1  and  a  model  atmosphere  from  Plez  (1992).  The  value  of  log  e  (N)  derived  by 
this  technique  is  not  a  true  nitrogen  abundance,  because  a  combination  of  C  and  N  abundances  are  needed  to  derive  a  true  value  of  log  e  (N)  from  CN  lines. 
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Fig.  6. — Relation  between  macroturbulent  velocity  ({)  and  surface 
gravity  (log  g)  for  supergiant  stars  (filled  circles).  All  supergiant  stars  from 
Table  1  and  IRS  7,  VR  5-7,  IRS  19,  and  IRS  22  from  Table  3  are  plotted. 
The  best  unweighted  linear  fit  to  the  data  ( solid  line)  is  used  to  compute  J 
for  the  GC  supergiant  stars  BSD  72,  BSD  124,  and  BSD  129. 


tainty  in  the  surface  gravity.  Equation  (3)  was  applied  to 
BSD  72,  BSD  124,  and  BSD  129. 

For  giant  stars,  the  range  in  log  g  is  too  narrow  to  derive 
a  relationship.  A  mean  value  of  £  =  9.0  ±  3.0  km  s“ 1  is  used 
for  the  GC  giant  stars.  This  mean  value  is  derived  from 
values  of  £  for  solar  neighborhood  giant  stars  and  IRS  11 
measured  by  the  broadening  of  the  CN  lines,  as  given  in 
Tables  1  and  3.  This  mean  value  was  adopted  for  BSD  114 
and  BSD  140. 

3.3.  Uncertainties 

Our  basic  observable  parameters,  Tat  and  Mbol,  have 
uncertainties  that  come  from  the  technique  used  to  compute 
them  (§§  3.2.1  and  3.2.2).  The  surface  gravity  uncertainty 
includes  both  the  uncertainty  in  Tat  and  the  uncertainty  in 
Mbol  (§  3.2.2).  The  microturbulent  velocity  uncertainty 
includes  both  the  uncertainties  in  log  g  and  the  uncertainty 
of  the  lit  used  to  calculate  it  (§  3.2.3).  The  same  is  true  for  the 
macroturbulent  velocity  (§  3.2.4). 

It  is  important  to  know  what  are  the  effects  of  each  of 
these  uncertainties  in  the  abundance  determination.  The 
uncertainty  in  [Fe/H]  for  each  stellar  parameter  has  been 
estimated  by  varying  one  of  the  stellar  parameters  and  com¬ 
puting  the  difference  in  iron  abundance.  A  typical  uncer¬ 
tainty  in  Teff  of  ±  300  K  results  in  an  uncertainty  of  +  0.08 
dex  in  [Fe/H],  A  typical  uncertainty  of  ±0.35  in  log  g 
implies  an  uncertainty  of  ±0.12  dex  in  [Fe/H].  A  typical 
uncertainty  of  ±0.4  km  s-1  in  £  causes  an  uncertainty  of 
+  0.09  dex  in  [Fe/H],  A  typical  uncertainty  of  ±  1.8  km  s-1 
in  £  translates  to  an  uncertainty  of  ±0.13  dex  in  [Fe/H], 
The  iron  abundance  is  not  very  sensitive  to  uncertainties  in 
Tat  and  £,  but  it  is  more  sensitive  to  uncertainties  in  log  g 
and  £.  In  addition  to  the  uncertainty  from  the  stellar  param¬ 
eters,  the  standard  error  must  also  be  considered  for  each 
star.  The  standard  error  comes  from  the  scatter  in  [Fe/H] 


as  derived  from  individual  Fe  i  lines.  All  the  uncertainties 
from  stellar  parameters  and  the  standard  error  for  each  star 
are  listed  in  Tables  7  and  8.  The  total  uncertainty  is  the 
quadratic  addition  of  the  uncertainties  derived  from  the 
stellar  parameters  and  the  standard  error  in  the  value  of 
[Fe/H]  from  individual  Fe  i  lines.  For  the  solar  neighbor¬ 
hood  stars  and  IRS  7,  the  quadratic  addition  of  the  uncer¬ 
tainties  is  a  good  estimation  of  the  total  uncertainty, 
because  all  the  uncertainties  come  from  independent  mea¬ 
surements.  In  the  cases  in  which  £  and  £  are  computed 
through  their  relationships  with  log  g,  the  uncertainties  are 
not  independent,  and  the  quadratic  addition  can  be  an 
overestimation  of  the  total  uncertainty.  We  estimate  that 
the  total  uncertainty  will  decrease  by  no  more  than  0.05 
dex,  if  we  consider  the  correlation  of  the  uncertainties  in 
log  g,  £,  and  £. 

Systematic  errors  may  be  present  as  well.  They  may  come 
from  unidentified  lines,  which  could  make  a  contribution 
that  is  not  included  in  our  synthetic  spectrum  because  of  the 
lack  of  atomic  and  molecular  parameters.  They  may  also 
come  from  errors  in  the  determination  of  (//-values  or  from 
failure  of  the  model  atmospheres  to  model  the  stars  we 
observe  correctly,  because  model  atmospheres  use  solar 
abundance  ratios,  or  NLTE  effects. 

3.4.  NLTE  Effects 

The  iron  abundance  could  be  affected  by  departures  from 
LTE.  The  main  NLTE  effect  in  late-type  stars  is  caused  by 
overionization  of  electron  donor  metals  by  ultraviolet  radi¬ 
ation  (Auman  &  Woodrow  1975).  Using  the  Saha  equation, 
we  found  that  in  our  typical  star  93%  of  iron  is  neutral. 
Thus  overionization  should  be  smaller  in  our  stars  than  in 
warmer  giants,  where  only  a  few  percent  of  iron  is  neutral. 

It  is  known  from  empirical  studies  that  abundances 
derived  from  low-excitation  lines  give  systematically  lower 
abundances  (~0.3  dex)  than  the  ones  derived  from  high- 
excitation  and  ionized  lines  in  late-type  stars;  high- 
excitation  and  ionized  lines  give  very  similar  abundance 
results  (Ruland  et  al.  1980;  Steenbock  1985;  Takeda  1991; 
Tomkin  &  Lambert  1983).  Recently,  Thevenin  &  Idiart 
(1999)  studied  NLTE  effects  in  Fe  abundances  in  metal- 
poor  late-type  stars.  They  found  that  ionized  lines  are  not 
significantly  affected  by  NLTE  and  that  NLTE  corrections 
become  less  important  as  [Fe/H]  increases,  being  minimal 
for  solar  abundance  stars.  Mostly  high-excitation  lines  are 
used  in  our  analysis,  so  the  NLTE  effects  that  might  be 
present  should  be  minimal. 

Tomkin  &  Lambert  (1983)  found  that  NLTE  effects 
cancel  out  when  a  differential  analysis  is  carried  out  relative 
to  a  very  similar  star  in  terms  of  effective  temperature  and 
luminosity.  Nearby  stars  of  similar  effective  temperature 
and  luminosity  as  the  GC  stars  (see  Fig.  2)  have  been 
analyzed  similarly  in  this  paper,  providing  a  differential 
abundance  comparison  that  should  remove  any  NLTE 
effects  that  might  be  present. 

4.  RESULTS 

Each  iron  abundance  value  for  a  particular  Fe  i  line  and 
star  was  derived  from  a  comparison  between  the  synthetic 
spectrum  generated  by  MOOG  and  the  observed  spectrum. 
The  [Fe/H]  results  for  each  line  for  cool,  luminous  stars  in 
the  solar  neighborhood  and  GC  stars  are  listed  in  Table  5 
and  6,  respectively.  The  uncertainties  per  line  are  estimated 
from  the  ability  to  distinguish  models  with  different  [Fe/H] 
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TABLE  5 

[Fe/H]  for  Nearby  Late-Type  Stars2 


Wavelength  of  Fe  i  Lines  (A) 


Star 

21781.8 

22381.3 

22386.9 

22391.2 

22399.0 

22818.8 

22838.6 

22852.2 

HR  6146  . 

-0.05 

+  0.05 

-0.20 

0.00 

-0.30 

+  0.15 

+  0.25 

+  0.10 

HR  6702  . 

-0.05 

-0.05 

-0.25 

-0.05 

-0.30 

+  0.20 

+  0.10 

+  0.30 

HR  7442  . 

+  0.10 

+  0.10 

-0.25 

+  0.10 

-0.60 

+  0.20 

+  0.20 

+  0.20 

HR  8062  . 

+  0.10 

-0.10 

-0.05 

+  0.05 

-0.20 

+  0.10 

-0.05 

a  Ori . 

+  0.15 

0.00 

+  0.05 

+  0.05 

-0.10 

+  0.20 

+  0.10 

0.00 

HD  202380 . 

+  0.10 

+  0.10 

+  0.10 

0.00 

+  0.10 

+  0.10 

+  0.10 

HR  8383  . 

-0.05 

-0.10 

-0.10 

-0.10 

-0.20 

-0.15 

-0.25 

-0.10 

HD  163428 . 

+  0.07 

0.00 

0.00 

-0.02 

-0.25 

-0.05 

-0.05 

BD  +59°594 . 

-0.10 

0.00 

-0.25 

-0.25 

-0.25 

-0.30 

-0.35 

-0.35 

HD  232766 . 

-0.10 

-0.20 

-0.15 

-0.15 

-0.40 

-0.20 

-0.40 

-0.40 

HR  8726  . 

+  0.15 

-0.10 

0.00 

-0.05 

-0.20 

0.00 

-0.10 

-0.20 

a  Error  in  [Fe/F[]  for  each  line  is  +0.05  dex;  [Fe/H]  determined  from  each  Fe  I  line. 


considering  the  S/N  of  the  observed  spectrum.  For  cool, 
luminous  stars  in  the  solar  neighborhood,  the  uncertainty 
per  line  is  estimated  to  be  ±  0.05  dex,  since  the  S/N  is  very 
homogeneous  among  all  lines  and  stars.  For  GC  stars,  the 
uncertainty  is  listed  individually  for  each  line  in  Table  6. 
Figures  1,  7,  and  8  show  the  synthetic  and  observed  spectra 
for  each  line  for  a  Ori,  cool,  luminous  stars  in  the  solar 
neighborhood,  and  GC  stars,  respectively. 

The  final  [Fe/H]  results  for  each  star  are  listed  in  Table  9. 
These  values  are  the  mean  [Fe/H]  weighted  by  the  individ¬ 
ual  uncertainties  in  [Fe/H]  for  each  Fe  i  line,  as  listed  in 
Tables  5  and  6.  The  final  uncertainty  in  Table  9  is  the  total 
uncertainty  from  Table  7  and  8. 

The  analyses  of  our  paper  and  Carr  et  al.  (2000)  have  five 
stars  in  common:  IRS  7,  a  Ori,  HR  6146  (30  Her),  HR  6702, 
and  HR  8383  (VV  Cep).  There  is  agreement  in  the  Fe  abun¬ 
dances  within  0.1  dex  for  four  stars  (IRS  7,  a  Ori,  HR  6146, 
and  HR  8383).  Our  Fe  abundance  for  HR  6702  is  0.18  dex 


lower  than  that  of  Carr  et  al.  (2000).  The  abundance  tech¬ 
nique  used  in  both  studies  is  very  similar,  but  Carr  et  al. 
(2000)  included  a  different  set  of  Fe  i  lines  and  slightly  differ¬ 
ent  stellar  parameters. 

IRS  7  is  known  to  be  a  variable  supergiant  (Blum  et  al. 
1996a).  There  is  one  Fe  line  common  to  the  analysis  of  Carr 
et  al.  (2000)  and  of  this  paper,  which  was  observed  at  differ¬ 
ent  epochs  by  Carr  et  al.  (2000)  and  by  us.  For  this  line,  our 
result  is  consistent  within  the  uncertainties  with  the  results 
of  Carr  et  al.  (2000).  Another  spectrum  of  IRS  7,  taken  with 
NIRSPEC  at  the  Keck  Telescope  1  year  after  our  obser¬ 
vations,  was  kindly  made  available  to  us  (D.  F.  Figer, 
private  communication).  We  have  one  Fe  line  in  common 
with  the  NIRSPEC  spectrum,  and  again  the  iron  abun¬ 
dance  results  from  both  data  sets  are  also  consistent  within 
the  uncertainties.  We  conclude  that  the  effect  of  variability 
on  our  IRS  7  results  are  negligible  considering  our  uncer¬ 
tainties. 


TABLE  6 

[Fe/H]  for  Galactic  Center  Stars* 


Wavelength  of  Fe  i  Lines  (A) 


Star 

21781.8 

22381.3 

22386.9 

22391.2 

22399.0 

22818.8 

22838.6 

22852.2 

IRS  7  . 

+  0.35 

+  0.22 

-0.15 

-0.10 

— 

+  0.38 

-0.30 

±0.05 

±0.07 

±0.10 

±0.05 

±0.05 

±0.05 

VR  5-7 . 

+  0.12 

+  0.20 

-0.30 

+  0.15 

-0.40 

+  0.10 

+  0.15 

±0.03 

±0.07 

±0.10 

±0.05 

±0.10 

±0.05 

±0.07 

IRS  19 . 

+  0.30 

+  0.40 

-0.30 

+  0.40 

-0.05 

+  0.45 

-0.05 

±0.05 

±0.05 

±0.10 

±0.07 

±0.10 

±0.05 

±0.10 

IRS  22 . 

+  0.30 

-0.05 

-0.40 

+  0.15 

-0.40 

+  0.20 

+  0.35 

+  0.30 

±0.05 

±0.05 

±0.10 

±0.10 

±0.10 

±0.10 

±0.10 

±0.10 

BSD  124 . 

+  0.22 

+  0.20 

+  0.20 

-0.10 

+  0.20 

+  0.20 

±0.07 

±0.10 

±0.20 

±0.10 

±0.05 

±0.10 

BSD  129 . 

+  0.50 

+  0.50 

+  0.70 

+  0.20 

+  0.60 

+  0.60 

±0.20 

±0.10 

±0.20 

±0.10 

±0.10 

±0.10 

BSD  72 . 

+  0.10 

+  0.40 

+  0.35 

+  0.05 

+  0.35 

-0.20 

±0.10 

±0.10 

±0.10 

±0.10 

±0.10 

±0.10 

BSD  114 . 

-0.20 

-0.85 

-0.10 

-0.70 

0.00 

-0.30 

±0.10 

±0.20 

±0.10 

±0.10 

±0.10 

±0.10 

IRS  11 . 

-0.15 

+  0.25 

-0.80 

-0.20 

-0.80 

-0.15 

-0.50 

±0.05 

±0.05 

±0.10 

±0.10 

±0.10 

±0.10 

±0.10 

BSD  140 . 

0.00 

-0.40 

+  0.10 

-0.40 

0.00 

0.00 

+  0.20 

+  0.20 

+  0.20 

+  0.20 

+  0.10 

+  0.10 

[Fe/H]  determined  from  each  Fe  i  line. 


Fe  21781.8  Fe  22381.3  Fe  22386.9  Fe  22391.2  Fe  22399  Fe  22818.8  Fe  22838.6  Fe  22852.2 


21780  22380  22386  22392  22398  22818  22836  22854 

x(A) 


Fe  21781.8  Fe  22381.3  Fe  22386.9  Fe  22391.2  Fe  22399  Fe  22818.8  Fe  22838.6  Fe  22852.2 


21780  22380  22386  22392  22398  22818  22836  22854 

x(A) 

Fig.  7. — Observed  (filled  squares)  and  synthetic  spectra  (solid  line )  for  cool,  luminous  stars  in  the  solar  neighborhood.  The  synthetic  spectrum  is  computed 
by  MOOG,  using  the  stellar  parameters  listed  in  Table  1  and  model  atmospheres  from  Plez  (1992).  The  error  bars  come  from  the  difference  of  observations 
taken  at  two  positions  along  the  slit.  Fe  lines  and  their  wavelengths  in  A  are  marked  (bold  vertical  lines )  at  the  top  panels.  CN,  Sc,  and  unidentified  lines 
(question  marks )  are  also  marked  (vertical  lines).  Tick  marks  along  the  x-axis  are  1  A  apart.  The  [Fe/H]  value  from  Table  5  for  each  line  is  given  in  each  panel. 
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TABLE  7 


Uncertainties  in  [Fe/H]  for  Solar  Neighborhood  Stars 


Star 

4-  T  a 

X  ieff 

±log  g“ 

±C 

±Ca 

Std.b 

Total' 

HR  6146  . 

+  0.03 

±0.10 

+  0.11 

±0.08 

±0.07 

±0.19 

HR  6702  . 

+  0.03 

±0.10 

+  0.11 

±0.13 

±0.08 

±0.22 

HR  7442  . 

+  0.03 

±0.10 

+  0.04 

±0.19 

±0.11 

±0.25 

HR  8062  . 

+  0.03 

±0.10 

+  0.04 

±0.04 

±0.04 

±0.13 

a  Ori . 

+  0.07 

±0.10 

+  0.04 

±0.04 

±0.03 

±0.14 

HD  202380 . 

+  0.05 

±0.17 

+  0.04 

±0.04 

±0.02 

±0.19 

HR  8383  . 

+  0.07 

±0.10 

+  0.04 

±0.15 

±0.03 

±0.20 

HD  163428 . 

+  0.05 

±0.17 

+  0.04 

±0.22 

±0.04 

±0.29 

BD  +59°594 . 

+  0.05 

±0.10 

+  0.04 

±0.07 

±0.05 

±0.15 

HD  232766 . 

+  0.05 

±0.10 

+  0.04 

±0.26 

±0.05 

±0.29 

HR  8726  . 

+  0.05 

±0.10 

+  0.04 

±0.04 

±0.04 

±0.13 

aTaf  =  effective  temperature;  g  =  surface  gravity;  {  =  microturbulent  velocity; 
f  =  macroturbulent  velocity. 

b  Standard  error,  determined  from  scatter  among  [Fe/H]  values  measured  for  different 
Fe  i  lines. 

c  Total  uncertainty,  derived  from  quadratic  sum  of  uncertainties  due  to  the  standard 
error  and  the  uncertainties  from  the  stellar  parameters  (Tat,  log  g,  (,  f). 


4.1.  Galactic  Center  Mean  [Fe/H] 

An  unweighted  mean  value  of  [Fe/H]  =  +0.09 
is  obtained  for  the  GC  stars.  A  mean  weighted  by  the  total 
uncertainties  in  [Fe/H]  (Table  9)  gives  a  value  of  [Fe/ 


H]  =  +0.12  for  the  GC  stars.  The  estimated  dispersion  or  a 
is  0.22  dex,  which  is  not  significantly  different  from  the 
typical  total  uncertainty  in  [Fe/H]  for  each  GC  star  of  0.28 
dex.  Even  if  the  total  uncertainties  are  overestimated  by 


TABLE  8 

Uncertainties  in  [Fe/H]  for  Galactic  Center  Stars 


Star 

+  T  a 

X  1  eff 

+  log  <7* 

±£a 

±1 

Std.b 

Total' 

IRS  7  . 

+  0.07 

±0.07 

+  0.09 

±0.19 

±0.13 

±0.27 

VR  5-7 . 

+  0.08 

±0.10 

+  0.11 

±0.11 

±0.11 

±0.23 

IRS  19 . 

+  0.08 

±0.10 

+  0.11 

±0.14 

±0.13 

±0.26 

IRS  22 . 

+  0.08 

±0.10 

+  0.11 

±0.11 

±0.12 

±0.24 

BSD  124 . 

+  0.08 

±0.10 

+  0.11 

±0.19 

±0.06 

±0.26 

BSD  129 . 

+  0.08 

±0.10 

+  0.11 

±0.19 

±0.08 

±0.27 

BSD  72 . 

+  0.08 

±0.10 

+  0.11 

±0.19 

±0.10 

±0.27 

BSD  114 . 

+  0.08 

±0.17 

+  0.13 

±0.21 

±0.16 

±0.35 

IRS  11 . 

+  0.08 

±0.17 

+  0.13 

±0.15 

±0.18 

±0.33 

BSD  140 . 

+  0.08 

±0.17 

+  0.13 

±0.21 

±0.10 

±0.32 

*  =  effective  temperature;  g  =  surface  gravity;  {  =  microturbulent  velocity; 

f  =  macroturbulent  velocity. 

b  Standard  error,  determined  from  scatter  among  [Fe/H]  values  measured  for 
different  Fe  i  lines. 

‘  Total  uncertainty,  derived  from  quadratic  sum  of  uncertainties  due  to  the  stan¬ 
dard  error  and  the  uncertainties  from  the  stellar  parameters  ( Ttff,  log  g,  {,  CI¬ 


TABLE  9 


Mean  [Fe/H]  for  Each  Star 


Solar  Neighborhood  Stars 

Galactic  Center  Stars 

Star 

^lines 

[Fe/H]s 

Star 

Clines 

[Fe/H]* 

HR  6146  . 

8 

-0.01  ±0.19 

IRS  7  . 

6 

+  0.09  +  0.27 

HR  6702  . 

8 

-0.02  ±0.22 

VR  5-7b . 

7 

+  0.09  +  0.23 

HR  7442  . 

8 

+  0.01  ±0.25 

IRS  19 . 

7 

+  0.29  +  0.26 

HR  8062  . 

7 

-0.03  +  0.13 

IRS  22 . 

8 

+  0.09  +  0.24 

a  Ori . 

8 

+  0.05  +  0.14 

BSD  124 . 

6 

+  0.17  +  0.26 

HD  202380 . 

7 

+  0.07  +  0.19 

BSD  129 . 

6 

+  0.49  +  0.27 

HR  8383  . 

8 

-0.14  +  0.20 

BSD  72 . 

6 

+  0.17  +  0.27 

HD  163428 . 

7 

-0.05  +  0.29 

BSD  114 . 

6 

-0.29  +  0.35 

BD  +59°594 . 

8 

-0.24  +  0.15 

IRS  11 . 

7 

-0.16  +  0.33 

HD  232766 . 

8 

-0.26  +  0.29 

BSD  140 . 

6 

-0.06  +  0.32 

HR  8726  . 

8 

-0.07  +  0.13 

*  Error  in  [Fe/H]  is  the  total  uncertainty  from  Tables  7  and  8. 
b  M  supergiant  star  in  the  Quintuplet  cluster  (R  =  30  pc). 
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~  0.05  dex  by  the  quadratic  addition  of  the  individual 
uncertainties,  the  typical  uncertainty  per  star  does  not  get 
small  enough  to  resolve  the  [Fe/H]  distribution  (see  §  3.3). 
This  means  that  the  dispersion  of  the  GC  [Fe/Fl]  distribu¬ 
tion  can  be  understood  by  the  uncertainties  present  in  the 
data. 

Note  that  the  supergiant  VR  5-7,  at  R  =  30  pc,  has 
[Fe/FI]  =  +0.09  ±  0.22,  which  is  similar  to  the  supergiants 
at  R  <  2.5  pc,  which  have  a  mean  [Fe/H]  =  +  0.21  ±0.15. 

An  unweighted  mean  of  [Fe/H]  =  —0.05  and  a  mean 
weighted  by  the  total  uncertainties  of  [Fe/H]  =  +0.03  was 
obtained  for  the  11  cool,  luminous  stars  in  the  solar  neigh¬ 
borhood.  The  estimated  dispersion  or  a  is  0.16  dex,  which 
again  is  not  significantly  different  from  the  typical  total 
uncertainty  in  [Fe/H]  which  is  0.20  dex  for  these  stars. 

The  obtained  mean  [Fe/H]  for  the  GC  is  close  to  solar 
and,  furthermore,  very  similar  to  the  cool,  luminous  stars  in 
the  solar  neighborhood  that  were  analyzed  in  the  same  way. 

5.  DISCUSSION 

Our  mean  [Fe/H]  for  the  GC  stars,  +0.12  +  0.22,  is 
similar  to  that  of  supergiants  in  the  Galactic  disk.  Figure  9 a 
compares  the  GC  [Fe/H]  distribution  with  the  [Fe/H]  dis¬ 
tribution  of  40  F,  G,  K,  and  M  supergiants  within  2  kpc  of 
the  Sun  (Luck  &  Bond  1989).  The  mean  [Fe/H]  of  the  solar 


-2  -1.6  -1.2  -0.8  -0.4  0  0.4  0.8  1.2 

[Fe/H] 


Fig.  9. — (a)  Fractional  distribution  of  [Fe/H]  for  10  GC  stars  ( solid 
line)  compared  to  [Fe/H]  for  40  solar  neighborhood  supergiant  stars 
(dashed  line)  from  Luck  &  Bond  (1989).  (b)  Fractional  distribution  of 
[Fe/H]  for  10  GC  stars  (solid  line)  compared  to  [Fe/H]  for  11  solar  neigh¬ 
borhood  stars  observed  and  analyzed  in  this  paper  (dashed  line),  (c)  Frac¬ 
tional  distribution  of  [Fe/H]  for  10  GC  stars  (solid  line)  compared  to 
[Fe/H]  for  262  Baade’s  window  stars  (dashed  line)  from  Sadler  et  al.  (1996). 


neighborhood  supergiant  stars  from  (Luck  &  Bond  1989)  is 
+  0.13  with  a  dispersion  of  0.20  dex.  The  width  and  mean  of 
the  two  distributions  agree  closely.  In  Figure  9b,  we 
compare  the  GC  [Fe/H]  distribution  with  that  of  the  11 
cool  supergiants  and  luminous  giants  we  have  analyzed  in 
this  paper.  Again,  the  mean  and  width  of  the  [Fe/H]  dis¬ 
tribution  for  the  GC  stars  and  solar  neighborhood  stars  are 
consistent.  Unlike  the  Luck  &  Bond  (1989)  sample  of  super¬ 
giants,  our  11  nearby  stars  are  restricted  to  the  same  range 
in  stellar  parameter  space  as  the  GC  sample ;  in  addition,  we 
have  used  the  same  analysis  method,  Fe  i  lines,  and  stellar 
models  for  all  the  stars.  It  is  important  to  emphasize  that 
this  differential  analysis  makes  the  abundance  comparison 
between  our  samples  of  GC  and  nearby  stars  robust  against 
possible  systematic  errors  that  might  effect  the  absolute 
values  of  [Fe/H],  Hence,  we  can  conclude  that  the  Fe  abun¬ 
dances  for  the  sample  of  GC  stars  in  this  paper  are  nearly 
identical  to  those  of  similar  Tef{  and  luminosity  stars  in  the 
solar  neighborhood,  with  no  evidence  for  a  supersolar 
metallicity. 

Our  mean  [Fe/H]  for  the  GC  stars  of  +0.12  +  0.22 
appears  to  be  in  conflict  with  the  iron  abundance  of  the 
Pistol  star,  at  R  =  30  pc  in  the  Quintuplet  cluster.  Najarro 
et  al.  (1999)  derive  [Fe/H]  ~  0.5  for  the  Pistol  star.  They 
qualify  their  result  as  preliminary,  since  details  of  the  atmo¬ 
sphere  modeling,  such  as  the  effects  of  charge  exchange 
reactions,  have  to  be  studied  and  included  in  future  work. 
The  Pistol  star  abundance  is  derived  from  emission  lines  in 
the  winds  from  this  hot  star,  so  the  atmospheric  modeling 
techniques  are  fundamentally  different  from  those  more 
established  techniques  we  have  used  to  study  abundances 
from  photospheric  absorption  lines  in  cool  giant  or  super¬ 
giant  stars. 

Serabyn  &  Morris  (1996)  have  proposed  that  the  central 
cluster  was  built  up  by  continuous,  perhaps  episodic,  star 
formation  over  the  lifetime  of  the  Galaxy.  In  this  case,  the 
central  100  pc  would  form  a  stellar  population  distinct  from 
the  old  population  of  the  Bulge.  In  Figure  9c,  the  GC 
[Fe/H]  distribution  is  compared  with  the  [Fe/H]  distribu¬ 
tion  for  262  stars  in  Baade’s  window  (BW)  in  the  bulge 
(l,  b  =  1°,  —4°;  Sadler,  Rich,  &  Terndrup  1996).  Sadler  et  al. 
(1996)  have  a  typical  uncertainty  per  star  of  +0.2  dex,  so 
their  finding  that  [Fe/H]  ranges  from  — 1.0  to  +0.5  in  13 W 
shows  the  abundance  spread  of  BW’s  distribution  is  well 
resolved  by  their  technique.  The  distribution  of  [Fe/H]  in 
the  GC,  however,  has  a  similar  width  as  the  uncertainties  in 
[Fe/H],  0.28  dex  per  star.  Stars  as  metal-rich  or  as  metal- 
poor  as  observed  in  the  Bulge  are  not  found  in  our  sample 
of  GC  stars.  However,  our  GC  stars  are  restricted  to  a 
young  to  intermediate-age  population,  and  hence  an  old 
stellar  population,  if  it  exists  in  the  central  100  pc,  is  not 
sampled  by  the  work  of  this  paper.  A  larger  abundance 
spread  could  be  represented  in  such  an  older  population, 
but  a  spectroscopic  abundance  analysis  will  require  the 
sensitivity  to  reach  K  and  M  giants.  A  more  complete  com¬ 
parison  must  also  include  measurements  of  the  abundance 
patterns  in  the  central  regions  relative  to  the  Bulge.  For 
example,  the  Bulge  K  giants  analyzed  by  McWilliam  & 
Rich  (1994)  show  a  distinct  enhancement  in  the  a  elements, 
compared  to  Fe. 

An  IMF  weighted  toward  more  massive  stars  in  the  GC 
has  been  proposed  by  Morris  &  Serabyn  (1996).  Recently, 
Figer  et  al.  (1999)  have  found  an  IMF  slanted  towards 
massive  stars  for  the  Arches  and  Quintuplet  clusters,  both 
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located  within  30  pc  of  the  GC.  A  history  of  chemical  evolu¬ 
tion  dominated  by  massive  stars  is  expected  to  result  in 
enhancements  of  oxygen  and  a  elements,  such  as  Mg,  Si,  Ca, 
and  Ti,  relative  to  Fe  (Wheeler,  Sneden,  &  Truran  1989). 
The  next  step  of  our  GC  project  is  to  measure  abundances 
of  a  elements  and  study  selective  enrichment  in  the  GC. 

Chemical  evolution  models  try  to  explain  abundance  pat¬ 
terns  by  considering  the  relative  star  formation  rate,  the  gas 
infall  and  outflow  rates,  the  star  formation  history,  and  the 
abundance  of  the  gas  compared  to  the  stars  (see  Audouze  & 
Tinsley  1976;  Rana  1991).  Most  chemical  evolution  models 
that  reproduce  the  radial  gradients  in  galaxies  are  valid 
only  for  distances  greater  than  2  kpc  from  the  Galactic 
center  (Tinsley  &  Larson  1978;  Samland,  Hensler,  &  Theis 
1997;  Chiappini,  Matteucci,  &  Padoan  2000;  Portinari  & 
Chiosi  1999).  Now  that  the  molecular  gas  in  the  GC  (from 
which  stars  form)  has  been  extensively  studied  (Stark  et  al. 
1991 ;  Blitz  et  al.  1993 ;  Morris  &  Serabyn  1996)  and  that  the 
abundances  of  stars  in  the  GC  has  been  obtained  for  the 
first  time,  we  strongly  urge  theoreticians  to  put  all  these 
data  together  into  a  detailed  chemical  evolution  model 
describing  the  central  parsecs  of  our  Galaxy. 

6.  CONCLUSIONS 

We  present  the  first  measurement  of  stellar  [Fe/H]  for  10 
stars  in  the  Galactic  center  (GC),  at  distances  from  the  GC 
of  R  <  30  pc.  Nine  GC  stars  are  located  in  the  central 
cluster  (R  <  2.5  pc),  and  one  star  is  located  in  the  Quin¬ 
tuplet  cluster  ( R  =  30  pc).  The  abundance  analysis  is  based 
on  high-resolution  (2/A2  =  40,000)  A -band  spectra.  The 
mean  [Fe/H]  of  the  GC  is  determined  to  be  near  solar, 
[Fe/H]  = +0.12  +  0.22,  and  also  similar  to  the  mean 
[Fe/H]  for  cool,  luminous  stars  in  the  solar  neighborhood, 
[Fe/H]  =  +0.03  +  0.16,  observed  and  analyzed  in  the 


same  way.  The  width  of  the  GC  [Fe/H]  distribution,  which 
ranges  from  [Fe/H]  =  —0.3  to  +0.5,  is  found  to  be  signifi¬ 
cantly  narrower  than  the  width  of  the  [Fe/H]  distribution 
of  Baade’s  window,  which  ranges  from  [Fe/H]  =  — 1.0  to 
+  0.8.  The  GC  [Fe/H]  distribution  is  consistent  with  the 
[Fe/H]  distribution  of  supergiant  stars  in  the  solar  neigh¬ 
borhood.  This  suggests  that  the  most  luminous  stars  in  the 
GC  are  unlikely  to  be  dominated  by  bulgelike  stars  and  that 
the  evolutionary  path  of  the  GC,  while  unique,  is  closer  to 
the  disk’s  than  to  the  bulge’s.  The  Quintuplet  star  at  R  =  30 
pc  has  a  similar  [Fe/H]  to  stars  located  in  the  central 
cluster  at  R  <  0.5  pc.  In  the  future,  abundance  measure¬ 
ments  of  CNO  and  a  elements  are  planned  to  provide  a 
complete  view  of  the  abundance  patterns  of  the  stars  in  the 
central  regions  of  the  Milky  Way. 
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